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We investigate a quasi one-dimensional ®^Rb Bose-Einstein condensate in a harmonic trap with an additional 
dimple trap (dT) in the center. Within a zero-temperature Gross-Pitaevskii mean-field description we provide 
a one-dimensional physical intuitive model, which we solve by both a time-independent variational approach 
and numerical calculations. With this we obtain at first equilibrium results for the emerging condensate wave 
function which reveal that a dimple trap potential induces a bump or a dip in case of a red- or a blue-detuned 
Gaussian laser beam, respectively. Afterwards, we investigate how this dT induced bump/dip-imprint upon 
the condensate wave function evolves for two quench scenarios. At first we consider the generic case that the 
harmonic confinement is released. During the resulting time-of-flight expansion it turns out that the dT induced 
bump in the condensate wave function remains present, whereas the dip starts decaying after a characteristic 
time scale which decreases with increasing blue-detuned dT depth. Secondly, once the red- or blue-detuned dT 
is switched off, we find that bright shock-waves or gray/dark bi-soliton trains emerge which oscillate within the 
harmonic confinement with a characteristic frequency. 

PACS numbers: 05.30.Jp, 32.80.Pj, 03.75.Lm, 05.45.Yv, 


I. INTRODUCTION 


The ability to manipulate and trap atoms with laser light has had a tremendous development in many fields of physics. The 
very first experimental success of trapping 500 Sodium atoms for several seconds in the tight focus of a Gaussian red-detuned 
laser beam occurred in 1986 JU]. The physical mechanism behind such an optical dipole trap is the electric dipole interaction 
of the trapped polarized atoms with the intense laser light, which is far detuned from the nearest optical transition of the atoms. 
They are hence largely independent from magnetic sublevels of the confined atoms, in contrast to a magneto-optical trap (MOT) 
which can only trap atoms with a certain internal state iBi. The so called dimple trap (dT) is nothing but a small tight 
optical dipole tr^ 141-1^. Cooling and trapping of atoms with these dT’s has a strong impact on the study of the Bose-Einstein 
condensates 10, Si. the observation of long decay times for atoms in their ground state li, and the research of trapping other 
atomic species or molecules M- 

A straightforward method for realizing a dT is to rely on the potential created by a freely propagating laser beam. The detuning 
of the laser frequency versus the atomic resonances determines, whether the atoms are red/blue-detuned, i.e. the laser frequency 
is below/above the atomic resonance frequency, respectively J^]. The red-detuned dT was particularly used for realizing matter 
wave traps EMI in the focus of a Gaussian laser beam. On the other hand, the blue-detuned Gaussian laser beam was used in 
optical waveguides EMI, where the creation of repulsive potentials was demonstrated^ using Laguerre-Gaussian laser beam 
iIMl. A focused or well-collimated Gaussian laser beam with a large red-detuning Eal or a dark hollow laser beam with a 
large blue-detuning 1271] were used to form 3D optical dipole trap’s, which can be widely applied to the accurate, non-contact 
manipulation and control of cold atoms fllMIli. 

In this paper, we will focus on studying neutral ®^Rb atoms within a quasi one-dimensional harmonic trap with an additional 
dimple trap. Experimentally, a highly elongated quasi-lD regime can be reached by tightly confining the atoms in the radial 
direction, effectively freezing-out the transverse dynamics O014381] . It is worth mentioning that, when the transverse length 
scales are of the order of or less than the atomic interaction len gth, the one-dimensional system can only be described within the 
Tonks-Girardeau or within the super-Tonks-Girardeau regime ]l39l - l4ll] . which is experimentally realizable near a confinement- 
induced resonance ll42l - l44l] . On the other hand, when the transverse confinement is larger than the atomic interaction strength, 
the underlying three-dimensional Gross-Pitaevskii equation (GPE) can be reduced to an effective quasi ID model ]|45[] . In one 
spatial dimension (ID) this equation is well-known, for instance, to feature bright and dark solitons for attractive and repulsive 
s-wave scattering lengths, respectively EMi. Many experiments investigate the collision of two Bose-Einstein condensates 
where the celebrated matter-wave interference pattern appears ]f^ or shock-waves are generated |5J.]. Eor lower collisional 
energies, the repulsive interaction energy becomes significant, and the interference pattern evolves into an array of gray solitons 
] 5214591] . Eurthermore, dark solitons can be created by manipulating the condensate density using external potentials 16014631] . 

This work is organized as follows; In Sec. II, we start with a model which describes the dynamical evolution of a quasi-lD 
Bose-Einstein condensate (BEC) in a magneto-optical trap with an additional red/blue-detuned dimple trap in the center. After- 
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wards in Sec. Ill, we justify a Thomas-Fermi approximation for the condensate wave function and compare it with numerical 
results. With this we show that the dT induces a bump or a dip upon the condensate wave function depending on whether dT 
laser beam is red- or blue-detuned. Subsequently, in Sec. IV, we discuss the dynamics of the dT induced bump/dip-imprint upon 
the condensate wave function for two quench scenarios. After having released the trap, the resulting time-of-flight expansion 
shows that the dT induced imprint remains conserved for a red-detuned dT but decreases for a blue-detuned dT. Furthermore, 
when the initial red/blue-detuned dT is switched off, we observe the emergence of bright shock-waves or gray/dark bi-soliton 
trains. Finally, Sec. V summarizes our findings for the proposed quasi-lD harmonically confined EEC with an additional dimple 
trap in the center in view of a possible experimental realization. 


II. QUASI ID MODEL 


We start with the fact that the underlying Gross-Pitaevskii equation for a condensate wave function can be formulated as the 
Hamilton principle of least action with the action functional 


.^30 = 



Esd (fir, 


where the Lagrangian density reads for three spatial dimensions 


^3D = 


ih 




dt 


dt 


2mB 


■0* (r, t) Alp (r, t)-V (r, t) ip (r, t) 




- Uij (r) iP* (r, t) Ip (r, f) - ^ || 0 (r, t) || 


( 1 ) 


Here ip (r, t) describes the EEC wave function with the coordinates r = {x, y, z) and the two-particle interaction strength reads 
Gg° = A^B47r/i^aB/mB, where Nb denotes the number of bosonic atoms. In case of ®^Rb atoms, the s-wave scattering length is 
Ob = 94.7 Rq with the Eohr radius Rq. We assume that the bosons are confined by a harmonic potential with an additional dT 
potential in the center. 

Eor instance, a MOT yields a harmonic confinement V (r) = mBUJ^z^/2 + tobWj^ -f y^) /2, which has rotational 
symmetry with respect to the z-axis. In the following, we consider the experimentally realistic trap frequencies Wr = 
27r X 160 Hz ^ Wz = Stt x 6.8 Hz l^l, so we have a cigar-shaped condensate, where the oscillator lengths amount to the 
values 4 = 0.84 pm <C 0 = 4.12 pm. 

An additional three-dimensional narrow Gaussian laser beam polarizes the neutral atoms and, thus, yields the dT potential 
[7® = UqI (r). Within the rotating-wave approximation its amplitude Uq is given by O6446711 


Sttc^F 




2u;lA 


( 2 ) 


Here F = | < e|d|p >1^ ui^/ (Sireofifi) denotes the damping rate due to energy loss via radiation, which is determined by the 
dipole matrix element between ground g and excited state e. Furthermore, A = oj — uja represents the detuning of the laser, 
where oj is the laser frequency and wa stands for the frequency of the Di- or D 2 -line of ®^Rb atoms, which are the transitions 
5^S'i/2 — 5^Pi/2 or 5^S'i/2 —>■ with the wave lengths Adi = 794.76 nm and Ad2 = 780.03 nm, respectively. The 

intensity profile of the Gaussian laser beam, which is assumed to move in y-direction, is determined via 


I(r) 


2P 

^Wx(y)VFz(y)® 




(3) 


where P denotes its power. Furthermore, Wx/z(y) = kkox/zy/l + ifi IVr^iz defines the beam radius in x- and z-direction, where 
the intensity decreases to 1/e^ of its peak value, and the Rayleigh len gths yRx/z = ttWox/z/A with wave length A = 27rc/a; 
define the distances where the beam radius increases by a factor of v/2 16511 . We use for the Gaussian laser beam width along 
the cc-axis = 1.1 pm and along the z-axis Woz = 3.2 pm, which are about ten times smaller than the corresponding ones 
used in Ref. Ij]. The corresponding Rayleigh lengths for the red-detuned laser light with A = 840 nm |[3l yield yRx = 4.526 pm 
and yRz = 38.29 pm and for the blue-detuned laser light with A = 772 nm we get yRx = 4.92 pm and yRz = 41.6 pm. Due 
to the fact yRx/z ^ h, we can approximate the widths of the beam in x- and z-direction according to Wx/z(y) ~ Wqx/z- This 
simplifies the dimple trap to 
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As the MOT provides a quasi one-dimensional setting due to ob h- Iz, we can follow Ref. Il45l] . and decompose the EEC 
wave-function ■!/)(r, f) = f)(^(rj_, f) with rj^ = {x, y) and 




97^ 






Subsequently, we integrate out the two transversal dimensions of the three-dimensional Lagrangian according to 


-ClD = 


/7 

J —OO j —GO 


£30 dxdy. 


(5) 


(6) 


After a straight-forward calculation the resulting quasi one-dimensional Lagrangian reads 


£id = 


ih 
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P;*iz,t)^^i^-Viz)tP*iz,t)tPiz,t) 


2m 
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(7) 


where V (z) = mBUJ^z^ j2 represents an effective one-dimensional harmonic potential from the MOT, and the one-dimensional 
two-particle interaction strength is 


Gb — ■ 

Furthermore, the one-dimensional dT depth turns out to be 

2UoP 


U = 


(8) 


(9) 


ttWoz^/WI + 21^ 

which depends on the power of the laser beam P as well as via (Ell and the detuning A = a; — wa on the laser wave length A. 
Note that the latter not only changes the absolute value of the dT depth but also its sign via the detuning A. For red detuning, 
i.e. when the laser frequency is smaller than the atomic frequency, the dT is negative and atoms are sucked into the dT potential. 
In the opposite case of blue detuning the atoms are repelled from the dT potential. Thus, the dT induces an imprint on the EEC 
wave function, which can be either a bump for red detuning or a dip for blue detuning. In the following we will analyze this 
interesting effect in more detail. 

To this end we consider the ID action 


Alin — 
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( 10 ) 


and determine the time dependent one-dimensional Gross-Pitaevskii equation (IDGPE) according to the Euler-Lagrangian equa 
tion 

SAw['ip*,'ip] dCiB d dCiD d dCm 


dz 

dz 


5tp* {z,t) dtp* [zA) 

Ey using the one-dimensional Lagrangian density (|7]i the IDGPE reads 


dt 

^ dt 
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d'^ 

2mB dz^ 


mjiUj- 2 

—-—z 


-fUe ^-f Gb II ||H0(z,f). 


( 11 ) 


( 12 ) 


On the right-hand side the first term represents the kinetic energy of the atoms with mass mB, the second term describes the 
harmonic MOT potential, the third term stands for the dT potential, and the last term represents the two-particle interaction. In 
order to make Eq. (fT2l i dimensionless, we introduce the dimensionless time as f = ujzt, the dimensionless coordinate z = zjlz, 
and the dimensionless wave function tp = tp/y/%. With this Eq. (fT2]) can be written in the form 


d 


1 52 


(57 = \ + T^z'^ + 'Oe -f Gb II 0 {z,i) f \tp{z,r), 


2dP 


(13) 


where we have Gb = 2NBU>raB/ (wz0) and U = U/ (huiz). For the above mentioned experimental parameters and Nb = 
20 X 10^ atoms of ®^Rb, we obtain the dimensionless couplings constant Gb = 11435.9. Furthermore, the typical dT depth 

|U| jkB = 210 nK yields the dimensionless value U = 643.83, and d = Wqz/ {V^lz) = 0.548 represents the ratio of the 
width of the dT potential along the z-axis and the longitudinal harmonic oscillator length. From here on, we will drop all tildes 
for simplicity. 
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FIG. 1: (Color online) Numerical density profile of BEC for the experimental coupling constanl value Gb = 11435.9 and for the dT depth U 
which increases from top to bottom according to the inlets. For a) negative values of U, the bump in the condensate wave function decreases, 
whereas for b) positive values the corresponding dip increases. 


X lO'' 



FIG. 2: (Color online) Ratio Sint,pot/F^kin versus U from solving IDGPE jl3t . 


III. DT INDUCED BUMP/DIP-IMPRINT UPON STATIONARY CONDENSATE WAVE EUNCTION 


In order to determine the dT induced imprint on the condensate wave function in equilibrium, we solve the IDGPE (fOl l 
in imaginary time numerically by using the split-operator method In this way we find that the dT-imprint leads to a 

bump/hole in the BEC density at the trap center for negative/positive values of U as shown in Eig. [T] Eor stronger red-detuned 
dT depth values the bump increases further, but for stronger blue-detuned dT the dip in the BEC density gets deeper and deeper 
until no more BEC atoms remain in the trap center. After this qualitative overview on the numerical results, we now work out 
an analytic approach for describing this red/blue-detuned dT induced bump/dip on the BEC density in a more quantitative way. 
To this end we present two arguments why the seminal Thomas-Eermi (TE) approximation is also applicable in our context. 

At first we provide a rough estimate in the case of an absent dT, i.e. U= 0, so the BEC density is characterized by the 
TE profile t/; (z) = \/{fJ. — z^/2) /Gb 0 (p — -z^/2), where the Heaviside function 0 prevents the density to become negative. 
Thus, the Tomas-Eermi radius i/2p follows from the dimensionless chemical potential /r, which is determined by normalization 
to be p = i As the red/blue-detuned dT is supposed to be inserted at the trap center, we then calculate the 

dimensionless BEC coherence length ^ at the trap center. It is defined by comparing the kinetic energy 1 with the interaction 
energy in the trap center, which is given by fi. Eor the above mentioned experimental parameters this yields the dimensionless 
BEC coherence length f = 0.038, which is about 14.4 times smaller than the dT width a = 0.548. This indicates that the dT 
induced imprint upon the BEC wave-function occurs on a length scale which is much larger than its coherence length, so the TE 
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approximation seems to be reasonable even in the presence of the red/blue-detuned dT. 

In view of a more quantitative justification for the applicability of the Thomas-Fermi approximation. Fig. |2] presents the 
numerical result how the ratio of the sum of interaction and potential energy i?int,pot versus the kinetic energy E]^in of the 
condensate wave function changes with increasing or decreasing the red/blue-detuned dT depth U. The maximal value of this 
energy ratio occurs for U= 0 and amounts to 7.5 x 10"^, which is of the order of the number of particles. Furthermore, we read 
off that the inequality iSint.pot/^'kin ^ 1 holds within the whole region of interest for U, so the TF approximation is there, indeed, 
valid. 

Therefore, we investigate in the following the TF approximation in more detail for non-zero red/blue-detuned dT depth U. To 
this end we use for the condensate wave function the ansatz ip{z, t) = insert it into the IDGPE (fl^ . and neglect the 

kinetic energy term, which yields the density profile 


^ = V ^ ■ T ■ 0 (m - Y - Ue-^) . (14) 

In view of the normalization \\ 'ip (z) \\'^ dz = 1, which fixes the chemical potential /i, we have to determine the Thomas- 
Fermi radii i?xF from the condition that the condensate wave function vanishes: 


F = 


r?2 

rtxF 


-Ue” 


(15) 


As can be read off from Fig.[T]the number of solutions of Eq. ( fTSl ) changes for increasing dT depth at a critical value Uc, which 
follows from solving the implicit equation 


Uc = i(^') (GB + v^aUe)t. 


(16) 


This yields the result Uc « 339.5 for the experimental coupling constant Gb = 11435.9, which compares well with the value 
Uc ~ 342 determined from solving IDGPE ( fTSl l. In the case that U is smaller than Uc Eq. ( fTST i defines only the cloud radius 
i?xFi- But for the case U > Uc the dT drills a hole in the center of the ®^Rb condensate, so it fragments into two parts. Thus, 
we have then apart from the outer cloud radius i?xFi also an inner cloud radius i?xF 2 - With this the normalization condition 

p-Rtfi 
7i?TF 2 


2 /ii™ II IP dz = 1 yields 


M (i?TFl - .RtF 2) - g (-^XFl “ -^TF2) “ - 




(17) 


where Erf(?/) = e ^ dx denotes the error function. In case of U < Uc the inner cloud radius i?xF 2 vanishes and the cloud 

radius is approximated via i?xFi ~ due to Eq. (fTSI l as it is much larger than the dimple trap width a. Thus, the chemical 
potential is determined explicitly from 




2/3 


(Gb -b v^aU) 


2/3 


U<Uc 


(18) 


Provided that U > Uc, the inner cloud radius i?xF 2 has to be taken into account according to Pig. [T] and, due to the fact that 


^TF 2 U, we get from Eq. (fTsT i the approximation p Ri Ue 


', which reduces to 


i?XF2 ~ at/log 


(19) 


Thus, we conclude that i?xF 2 vanishes, indeed, at Uc according to Eq. (fTbl l and Eq. (fTsT l. With this we obtain from Eq. (fTTl) that 
the chemical potential follows from solving 


v^aU -b Gb -b 2afiJlog ^ 


3v^Q;UErf 


( V log (^) ) + ^og^ + 4^p3/2 ^ u > Uc. (20) 


Pigure [3] shows the resulting outer and inner Thomas-Permi radius as a function of the dT depth U. We read off that T/xfi ~ 
1/271 remains approximately constant for U > Uc, so we conclude that the chemical potential p is locked to its critical value 
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U 

FIG. 3: (Color online) Outer Thomas-Fermi radius 7?tfi (red solid) and inner Thomas-Fermi radius i?TF2 (blue dashed) versus dimple trap 
depth U. EEC fragments into two parts above Uc « 339.5. 



FIG. 4: (Color online) Condensate density for U = 1000 from solving IDGPE dl3b in imaginary time numerically (red) and from TE 
approximation d (blue-circles). 


/ic ~ Uc = 339.5. Furthermore, we note that the inner Thomas-Fermi radius i?TF 2 increases up to about 5.4 q; for the considered 
range of U. 

Figure|4]compares the resulting TF condensate wave function (fT4l i with a numerical solution of the IDGPE (fOl l in imaginary 
time at U= 1000 and we read off that both agree quite well. Thus, our TF approximation describes the equilibrium properties of 
the condensate wave function in the presence of the red/blue-detuned dT even quantitatively correct. In view of a more detailed 
comparison, we characterize the red/blue-detuned dT induced imprint upon the condensate wave function ip{z) by the following 
two quantities. The first one is the hight/depth (HD) of the dT induced imprint 

r iiV'Wiis-iiV'WiiLo u<o 

HD = (21) 

[ Max(|| 2 /)( 2 ;) ||2) - II ,/,(0) ||2 U>0 

and the second one is the red/blue-detuned dT induced imprint width W, which we define as follows. For U < 0 we use the full 
width half maximum 


II ^ (W/2) \\l= (II V' (0) 115 + II ^ (0) ||5^o) /2 U < 0, 
whereas for U > 0 we define the equivalent width llll]: 


' -^iviax \ 

2 /oZMax- J \\-tp {z) \\l dz\ / {lo- \\ Ip {0) \\^) U>0, 


( 22 ) 


W = 


(23) 
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FIG. 5: (Color online) a) Height/depth and b) width of the dT induced hump/dip according to Eqs. dzB-lll, versus the red/blue-detuned dT 
depth U for the experimental EEC coupling constant Gb = 11435.9 calculated numerically by solving IDGPE ( 113b in imaginary time (blue 
circles) and analytically (black stars) from the TE condensate wave function (O. 


where we have Iq = Max (|| ip {z) ||u). Figurej^a) shows the red/blue-detuned dT induced imprint height/depth as a function 
of U. At first, we read off that for U = 0, i.e. when we have not switched on the dT, bump/dip vanishes. Furthermore, in the 
range U < Uc we observe that height/depth of the dT induced imprint bump/dip changes linearly with U according to 

|U| 

HD « y . (24) 

Gb 

In case of U > Uc height/depth of the dT induced imprint has approximately the constant value HDc = Uc/Gb ~ 0.029 as 
follows from the TF wave function (fTdl i and the above mentioned locking of the chemical potential to its critical value. Note that 
this constant value only slightly deviates from the corresponding numerical value HDc = 0.03. 

Correspondingly, Figure |5]b) depicts the dimple trap induced width W as a function of U. From our TF approximation we 
obtain for the width transcendental formulas, which read in case of U < 0 


—+ 2Ue 4. _u+-(-j 


2/3 


- (Gb + v^aU) 


2/3 


= 0 , 


and forU > o? j 2 


W = 


2a^ Y^log (^) [2 log (^) + 3] — fiv^aUErf ^^log 


3[a2 + a2log(|y) -2U] 


(25) 


(26) 


As shown in Fig.|5]b), for an increasing red-detuned dT depth, the width remains approximately constant, but just before U = 0 
starts to decrease to zero. For a blue-detuned dT the width of the dip continuously increases with an intermediate plateau at Uc 
with the value Wc ~ 1.91, which agrees well with the numerically obtained one Wc ~ 1.99. 


IV. DIMPLE TRAP INDUCED IMPRINT UPON CONDENSATE DYNAMICS 

In an experiment, any dT induced bump/dip upon the condensate wave function could only be detected dynamically. Thus, it 
is of high interest to study theoretically whether the dT induced bump/dip imprint, which we have found and analyzed for the 
stationary case in the previous section, remains present also during the dynamical evolution of the condensate wave function. 
To this end we investigate two quench scenarios numerically in more detail. The first one is the standard time-of-flight (TOF) 
expansion after having switched off the harmonic trap when the amplitude of the red/blue-detuned dT is still present. In the 
second case we consider the inverted situation that the red/blue-detuned dT is suddenly switched off within a remaining harmonic 
confinement, which turns out to give rise to the emergence of bright shock-waves or bi-solitons trains, respectively. 
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FIG. 6: (Color online) Height/depth of the dT induced imprint after having released the harmonic trap versus time for a) increasing negative 
and h) decreasing positive values of dT depth U from top to hottom. Inlet; relaxation time fid decreases with increasing U. 


A. Time-of-Flight Expansion 

Time-of-flight (TOF) absorption pictures represent an important diagnostic tool to analyze dilute quantum gases since the 
field’s inception. By suddenly turning off the magnetic trap, the atom cloud expands with a dynamics which is determined by 
both the momentum distribution of the atoms at the instance, when the confining potential is switched off, and by inter-atomic 
interactions if^l^ . We have investigated the time-of-flight expansion dynamics of the BEC with the dT by solving numerically 
the IDGPE (fTSl l and analyzing the resulting evolution of the condensate wave function. It turns out that, despite the continuous 
broadening of the condensate density, its dT induced imprint remains qualitatively preserved both for red and blue-detuned dT. 
Therefore, we focus a more quantitative discussion upon the dynamics of the corresponding dT induced imprint height/depth 
and width. 

Eor a red-detuned dT, it turns out that the bump height even remains constant in time. This is shown explicitly in Eig. |6]a), 
which roughly preserves its initial value at f = 0. In case of the bump width, we even find that no significant changes do occur 
neither in time nor for varying U, therefore we do not present a corresponding figure. Note that the latter finding originates from 
Eig.|5]b), where the width is shown to be roughly constant for all dT depths. 

Instead, in case of a blue-detuned dT, the dip decays after a characteristic time scale as shown in Eig.|6]b). The inlet reveals 
that the dip relaxes with a shorter time scale for increasing blue-detuned dT depth U. In addition, we read off from Eig. |7] that, 
at the beginning of TOE, the dT induced imprint width remains at first constant and then increases gradually. This change of W 
occurs on the scale of the relaxation time of HD, which is depicted in Eig.|6]b). 


B. Wave Packets Versus Solitons 


Due to their quantum coherence, BECs exhibit rich and complex dynamic patterns, which range from the celebrated matter- 
wave interference of two colliding condensates M\ fsill over Earaday waves 1173 .It^I to the particle-like excitations of solitons 
1 52-^ ^ W]. Eor our ID model of a BEC with a harmonic and a dimple trap in the center, we investigated the dynamics of 
the condensate wave function which emerges after having switched off the dT. To this end, Eig. [8] depicts the resulting profile 
of density n = |'0p and phase (p = tan“^ (V’Re/V'im) of the condensate wave function ip at different instants of time. Both 
for an initial red- and blue-detuned dT, we observe that two excitations of the condensate are created at the dT position, which 
travel in opposite direction with the same center-of-mass speed, are reflected at the trap boundaries and then collide at the dT 
position. Eurthermore, we find that these excitations qualitatively preserve their shape despite the collision and that the BEC wave 
function reveals characteristic phase slips between —7r/2 and tt/2. All these findings are not yet conclusive to decide whether 
these excitations represent wave packets in the absence of dispersion or solitons. Therefore, we investigate their dynamics in 
more detail, by determining their center-of-mass motion via 


S-Zq ^ (II ^ llu ~ II ^ llu=o) 
f-Zo (II ^ (^’ llu “ II ^ (^’ llu=o) 


^^L.R (t) 


(27) 
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FIG. 7: (Color online) Width of the dT induced imprint after having released the trap versus time for decreasing positive values of dT depth U 
from top to bottom. 


which are plotted in Fig. |9] Note that the mean positions zl and zb. of the excitations are uncertain in the region where they 
collide. Nevertheless Fig.|9]demonstrates that the excitations oscillate with the frequency il = 2 tt x 4.87 Flz irrespective of sign 
and size of U. As we have assumed the trap frequency = 27r x 6.8 Hz, we obtain the ratio H/wz ~ 0.72, which is quite close 
to H/wz = 1/V2 « 0.707. 

Despite these similarities of the cases of an initial red and blue-detuned dT, we observe one significant difference. Whereas 
the oscillation amplitudes of the excitations do not depend on the value of the initial U < 0 according to Fig. |9] a), we find 
decreasing oscillation amplitudes of the excitations with increasing the initial U > 0 in Fig.|9]b). Such an amplitude dependence 
on the initial condition is characteristic for gray/dark solitons according to Ref. ifs^ . This particle-like interpretation of the 
excitations agrees with the other theoretical prediction of Ref. if^ that gray/dark solitons oscillate in a harmonic confinement 
with the frequency H/wz = l/v/2, which was already confirmed in the Hamburg experiment of Ref. S and is also seen in 
Fig.|9] 

Conversely, for an initial red-detuned dT the excitations can not be identified with bright solitons as the dynamics is governed 
by a GPE with a repulsive two-particle interaction. Here the excitations have to be interpreted as wave packets which move 
without any dispersion as follows from a Bogoliubov dispersion relation and the smallness of the coherence length. Thus, for 
U < 0 the excitations propagate like sound waves in the BEC and, within a TE approximation, their center-of-mass motion 
is described by the evolution equation 17711 


dz{t) I 
dt ~ y ^ ^ 


(28) 


Solving ( 1281 ) with the initial condition ^(0) = 0 yields the result z{t) = sin f/v/2. Thus, we read off that the oscillation 

amplitude coincides with the TE radius and that the dimensionless oscillation frequency turns out to be H = l/v/2 in agreement 
with Eig.|9]a). 

Thus, we conclude that switching off the red/blue-detuned dT leads to physically different situations. Eor an initial red-detuned 
dT, we generate wave packets which correspond to white shock waves 17811 . whereas for the corresponding blue-detuned case 
bi-soliton trains emerge 58, 79], due to the collision of the two fragmented parts of the condensate. Note that it can be shown 
in our system that gray bi-solitons trains are generated for a partially fragmented BEC, i.e. U < Uc. On the other hand the dark 
bi-solitons trains turn out to be only generated for U > Uc, where the BEC is well fragmented into two parts equilibrium. 


V. SUMMARY AND CONCLUSION 

In the present work we studied within a quasi ID model both analytically and numerically how a dimple trap in the center 
of a harmonically trapped BEC affects the condensate wave function. At first, we showed for the equilibrium properties of 















10 


t=0.45 t=0.75 t=l.l t=3.6 t=5.2 




FIG. 8: (Color online) Density (phase) profile of EEC after having switched off the red/blue-detuned dT: in blue-dashed line at t=0, and in 
red-solid line at f=0.45 (1'’* column), t=0J5 (2"'^ column), f=l.l (3'^'^ column), f=3.6 (4*^ column), and t=5.2 (5**' column) for (a) [/=-100 
and (b) (7=100. 




0 2 4 6 8 10 

t 


FIG. 9: (Color online) Center of mass positions of excitations Zh (filled circles) and zr (empty circles) according to Eq. versus time after 
having switched off the dT with increasing absolute value of the depth |U| from top to bottom, for a) red-detuning and b) blue-detuning. Black 
filled circles represent the region of colliding excitations, where mean positions are not perfectly detectable. 


the system that the Thomas-Fermi approximation agrees quantitatively with numerical solutions of the underlying ID Gross- 
Pitaevskii equation. For an increasing red-detuned dT depth, it turns out for the induced bump that the height decreases linearly, 
whereas the width remains approximately constant. In contrast to that we found for an increasing blue-detuned dT that depth 
and width of the induced dip initially increase. Beyond a critical value Uc, the BEC even fragments into two parts and, if U is 
increased beyond Uc, the dT induced imprint yields a condensate wave function whose width increases further, although the dip 
height/depth remains constant. Afterwards, we investigated the dT induced bump/dip upon the condensate dynamics for two 
quench scenarios. 
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At first, we considered the release of the harmonic confinement, which leads to a time-of-flight expansion and found that the 
dT induced imprint remains conserved for a red-detuned dT but decreases in the blue-detuned case. This result suggests that 
it might be experimentally easier to observe the bump for a red-detuned dT. On the other hand, in an experiment one has to 
take into account that inelastic collisions lead to two- and three-body losses of the condensate atoms IS^, 811- As such inelastic 
collisions are enhanced for a higher BEC density, they play a vital role for a red-detuned dT, when the condensate density has 
a bump at the dT position, but are negligible for the blue-detuned dT with the dip in the BEC wave-function. Thus, a more 
realistic description of the experiment needs to consider the loss of condensate atoms by adding damping terms to the IDGPE 
O, which are of the form 1 X 2 || 'k {z) |p and iT^ || {z) 11 "^, where the positive constants T 2 and T 3 denote two- and 

three-body loss rates, respectively. We note that these additional terms may have nontrivial effects on the dT properties 0 . 

In addition, we analyzed the condensate dynamics after having switching off the red/blue-detuned dT. This case turned out 
to be an interesting laboratory in order to study the physical similarities and differences of bright shock-waves and gray/dark 
bi-soliton trains, which emerge for an initial red- and blue-detuned dT, respectively. The astonishing observation, that the oscil¬ 
lation frequencies of both the bright shock-waves and the bi-soliton trains coincide, is presumably an artifact of the harmonic 
confinement. Thus, it might be rewarding to further investigate these different dynamical features also in anharmonic confine¬ 
ments Additionally, we have also found that the generation of gray/dark bi-soliton trains is a generic phenomenon on 

collisions of partially/fully fragmented BEC, respectively, and the partially/fully fragmented BEC is strongly depending upon 
the equilibrium values of the dimple trap depth. 
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